Ultracold Thermalization of 7 Li and 87 Rb 
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We report on measurements of cross-species thermalization inside a magnetically trapped spin- 
polarized mixture of 87 Rb and 7 Li atoms with both atoms in their respective low field seeking 
magnetic substates \F = 2,mF = 2). Measurement of the thermalization velocity in the ultracold 
regime below 10 /iK allows for the derivation of the absolute value of the pure triplet s-wave scattering 
length governing the interaction. We find |ar,sr| = (59 ± 19) o,b- We propose to study both species 
in the condensed regime to derive the sign of 0737. In this context, we present numerical solutions 
to the coupled Gross-Pitaevskii equation based on the hypothesis of a positive sign. According 
to the simulations, phase separation of the Li and Rb |2, 2) clouds occurs along with a mean- 
field stabilization allowing for larger atom numbers of condensed 7 Li atoms before collapse sets in. 
Observation of this mean-field stabilization would directly fix the sign of (17,87. We discuss our 
results in the light of this proposal. 

PACS numbers: 34.50.-s, 03.75.Mn, 32.80.Pj 
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I. INTRODUCTION 

Sympathetic cooling with an independent heat bath 
is a convenient tool to cool atomic clouds that cannot 
be cooled directly. Its application to Fermi-Bose mix- 
tures led to a breakthrough in reaching the quantum de- 
generate regime with fully spin polatized fcrmionic atom 
clouds. The antisymmetry of the fermionic wave func- 
tions prohibits direct evaporative cooling due to miss- 
ing s-wave scattering. Adding a second bosonic species 
that can have collisions with the fermionic species cir- 
cumvents this difficulty. A range of different Fermi-Bose 
mixtures was brought to the degenerate regime employ- 
ing this technique: 6 Li+ 7 Li 0], 6 Li+ 23 Na @, 40 K+ 87 Rb 
0, 6 Li+ 87 Rb 3, 3 Hc+ 4 He @. For bosonic target 
species it is also interesting to be cooled by an inde- 
pendent cooling agent. Two 87 Rb condensates in dif- 
ferent hyperfine states F could be produced by evapora- 
tion of atoms in the F = 2 state while atoms in F = 1 
were passively cooled through thermalizing collisions Q . 
Another example is 85 Rb, whose intra-species scatter- 
ing length ass has a zero crossing at collision energies 
Ecoii/ks ~ 375 /iK. This prohibits direct evaporation 
of 85 Rb starting with higher temperatures. The prob- 
lem was overcome by sympathetically cooling 85 Rb in a 
bath of 87 Rb p}. Sympathetic cooling down to quan- 
tum degeneracy of the target species was also realized 
in 41 K+ 87 Rb and recently in 39 K+ 87 Rb Q3 where 

again 87 Rb served as coolant. Furthermore, interspecies 
scattering properties were determined in a mixture of 
87 Rb and 133 Cs, which was brought to low temperatures 
with the same technique [ll|. Mixtures of different iso- 
topes of Yb Q3] as well as a 174 Yb+ 87 Rb mixture [3 
are currently under examination as well. 

When working with 7 Li, it is important to note that 
Bose-Einstein condensates of this species are intrinsically 
unstable in the \F, mj?) = |2,2) hyperfine state because 
the 7 Li triplet scattering length a 7 = — 23.7 <2b is nega- 



tive. This instability and the maximum condensate size 
were predicted theoretically [l4[ and later confirmed ex- 
perimentally [l5[. The strong effects that negative scat- 
tering lengths have on condensates were also shown in the 
controlled destabilization [l6[ and the spectacular 'bose- 
nova' experiments (l7| in Bose-Einstein condensates of 
85 Rb, where the strength and sign of the scattering length 
was rapidly changed by means of a Feshbach resonance. 
Our interest in cooling 7 Li with 87 Rb originates in the 
mean-field interaction in the regime where both species 
are Bose-condcnscd. Depending on the strength of this 
interaction, it can strongly alter the stability conditions 
for 7 Li Bose-Einstein condensates, as we show in a theo- 
retical section of this paper. 

Here, we report on the experimental realization of sym- 
pathetic cooling of 7 Li in a bath of 87 Rb. From controlled 
thermalization experiments in a temperature range from 
10 /iK to 2 /iK we extract the interspecies scattering cross 
section 07.87 = <iira% 87 . This fixes the absolute value of 
the triplet scattering length to 07,87! = (59 ± 19) ae • 
Upon cooling further down, we reach criticality with re- 
spect to Bose-Einstein condensation in 7 Li before 87 Rb 
starts to condense. The 7 Li criticality manifests itself in 
additional sudden losses of the 7 Li atom number as for 
our trap conditions the atom number of the overcritical 
7 Li largely exceeds the atom number of stable 7 Li con- 
densates. 

We also present numerical solutions of the coupled 
Gross-Pitaevskii equations of two interacting condensates 
of 7 Li and 87 Rb in spherically symmetric traps. Stabiliza- 
tion of the Li condensate for large Rb atom numbers is 
observed for a positive sign of 07,87- Experimentally, ob- 
servation of this effect can be expected when intra- and 
interspecies inelastic loss coefficients are small and the 
overlap is good. A precondition for this is that the rela- 
tive gravitational sag is smaller than the Thomas-Fermi- 
radii of the condensates involved. Hence, by observing 
the collapse of Li and the impact of the mean-field en- 
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ergy of Rb on the collapse in a regime where both species 
are condensed, it is, in principle, possible to determine 
the sign of the mutual interaction [18| ■ 

We structure the paper as follows: In Section [TT] we 
shortly describe the apparatus and the experimental pro- 
cedure used for cooling. In Section [ill! experimental re- 
sults for crossed species thcrmalization are presented and 
the absolute value of the scattering length determined. 
Section IIVI presents experimental evidence that we are 
able to reach the critical temperature for Bose-Einstein 
condensation of 7 Li. Finally, Section [V] describes the nu- 
merical calculations for coupled 7 Li and 87 Rb conden- 
sates. 



tical absorption imaging. The trap is rapidly turned off 
and the clouds expand ballistically for a variable amount 
of time. Light pulses resonant with the F = 2 — > F' = 3 
transitions of the respective £>2-hnes then map the den- 
sity distribution onto CCD cameras. In the case of Li, 
an additional counter-propagating rcpumping beam res- 
onant with the F = 1 — > F' = 2 transition is applied to 
keep the atoms in the cycling transition for the absorp- 
tion imaging. The purity of the samples with respect 
to Zeeman substates is monitored by application of a 
magnetic field gradient during ballistic expansion of the 
clouds. 



II. EXPERIMENTAL SETUP 

Our experimental setup is identical to the one de- 
scribed in [3| where we studied thermalization between 
fermionic 6 Li and 87 Rb, except for the fact that we now 
use the bosonic isotope of Li. As atom sources we use 
commercial dispensers with (enriched) 87 Rb and a Zee- 
man slower for the 7 Li. Up to 3 • 10 8 Rb and 8 ■ 10 7 Li 
atoms are simultaneously trapped in two spatially over- 
lapping magneto-optical traps (MOT). We then compress 
the MOTs by increasing the trap coil currents and then 
apply a short (1 ms) stage of polarization gradient cool- 
ing. The Rb then has a temperature of T&j ps 80 /uK. 
The Li, for which polarization gradient cooling does not 
work because of the unresolved hyperfine structure in 
the excited state, has an estimated temperature of sev- 
eral hundred [iK. Before turning on the magnetic trap, 
the Rb is optically pumped into the stretched state |2, 2) 
while Li is left unpolarizcd. We magnetically trap al- 
most all Rb atoms and approximately 10 7 Li atoms fl9j |. 
The mixture is then adiabatically transferred via several 
intermediate magnetic traps into a Ioffc-Pritchard type 
trap characterized by the secular frequencies ujsi/2tt — 
(50 x 206 x 200) 1 / 3 Hz for the Rb |2, 2) st ate (fo r Li |2, 2) 
the frequencies are higher by a factor y/87/7) and the 
magnetic field offset 3 G. As in the case of 6 Li, purity of 
the Rb spin states has to be ensured by special measures. 
Directly after turning on the Ioffe-Pritchard trap, a clean- 
ing microwave sweep coupling residual Rb |1, — 1) atoms 
to the untrapped |2, —2} state is applied. The Rb cloud is 
then selectively cooled by forced microwave evaporation 
coupling the |2, 2) state with the anti-trapped |1, 1) state 
at a certain distance from the trap centre. An additional 
microwave knive has to be applied that constantly emp- 
ties the trap from 2, 1) Rb atoms by coupling them to the 
untrapped |1,0) state. The |2, 1) atoms accumulate dur- 
ing the last stages of evaporation via dipolar relaxation 
[20j and lead to severe losses in Li atom numbers through 
inelastic collisions. The Li cloud, which is not actively 
cooled, adjusts its temperature to the Rb cloud via in- 
terspecies thcrmalization. We find that 7 Li thermalizes 
much faster than its fermionic counterpart Li Q. Atom 
numbers and temperatures are derived by standard op- 



III. CROSS SPECIES THERMALIZATION 

To measure the interspecies thermalization speed, we 
perform experiments at three different temperatures. 
Each experiment consists of two stages. First, the mix- 
ture is slowly (20 s) cooled to a common temperature 
T\. Then, a stage of rapid cooling (50 ms) follows, which 
leaves the Rb at a temperature T2, thus producing a tem- 
perature difference between the two atomic ensembles. 
The thermalization of the Li in the colder Rb bath is used 
for determining the scattering length using a straight for- 
ward collision model [2l[ . The temperature ranges in the 
three measurements are Ti//iK— ► T^j fiK: 10 — ♦ 6, 6 — > 4 
and 4.5 — ► 2. On a timescale of 300 ms thermalization 
is completed, which is considerably faster than for 6 Li, 
where thermalization times on the order of 3 s were ob- 
served under similar conditions [J]. 

It follows a short description of the thcrmalization 
model which we use to derive the scattering length. The 
total collision rate depends on the overlap integral of the 
two atomic clouds which, for purely thermal samples in 
harmonic traps, can be exactly solved [22| 
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with the collision cross section 07.87 = 4-770,7 g7 , the m " 
terspecies pure triplet scattering length 0,7,87, the mean 
trap frequency C07 = (uj^ r W7 iZ ) 1 / 3 and the mean rela- 
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with equal masses, an average of 2.7 collisions is needed 
for thermalization [23[, however, in the case of unequal 
masses efficiency is reduced by a factor £ = 4m 7 m g7 

j22| . In our case, Li needs ~10 collisions to thermal- 
ize with the Rb. Here we assume the Rb temperature 
not to change during the thcrmalization process. We en- 
sure this experimentally by application of a constant fre- 
quency microwave knife during the thermalization stage 
of the experiment which keeps the Rb temperature (but 
not the atom number) fixed. The temperature difference 
AT = T7 — Tg7 follows the differential equation 



d(AT) 
dt 
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with the collision rate per Li atom = . 
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FIG. 1: Measured (data points) and simulated evolution 
(solid curve) of the thermalization process in the case of 
21 = 10 /iK. After about 300 ms the thermalization process 
is complete. Each datapoint is an average of three individual 
measurements, the errorbars indicate the statistical uncer- 
tainty. 

Application of this model to the measured data results 
in values for |<ir ; 87| summarized in table |TJ 



Ti/fjK 


T 2 /fiK 


«7,87 |/OB 


x 2 


10 


6 


56 


0.500 


6 


4 


65 


0.238 


4.5 


2 


56 


0.924 



TABLE I: Resulting absolute scattering lengths for three dif- 
ferent thermalization experiments. The quality of the fit \ 2 
is also shown. The variations in the derived scattering length 
shows the limitation as to the accuracy in determining |ar,sr| 
with this method. The mean value is a — 59 sb with a stan- 
dard deviation es = 5.2 as- Further errors are discussed in 
the main text and Fig|2] 

A correction has to be applied to the ballistic expan- 
sion time for the Li component. This is because the self- 
inductance of the trap coils results in a coil current shut- 
down timescale on the order of 200 /is. The cloud of Li 
atoms already starts to expand during this turn off pro- 
cess whereby the effective time-of-flight is larger than the 
nominal value. Hence, the temperatures derived from 
the Gaussian widths of the Li clouds using the nomi- 
nal time-of-flight value are overestimated. Therefore, an 
additional experimentally determined correction to the 
time-of-flight of 130 /is is added to the nominal time of 
flight of 1000 /is. This changes the measured tempera- 
tures by up to 25 % for low temperatures as compared 
to uncorrected values. For very slow evaporation ramps, 
where Li and Rb are in constant equilibrium, this correc- 
tion yields the same values for both temperatures. The 
delayed turn off effect is negligible for the heavier Rb 
atoms which we let expand for a much longer time (20 
ms). 



The error introduced due to inexact values for the 
temperature difference can be estimated by applying a 
small temperature offset to the measured temperature 
evolution. The effect on the derived scattering length 
for Ti = 10 /iK is shown in Fig. [2k,. Mismatched tem- 
perature estimates of only 200 nK in either direction re- 
sult in variations in |<Z7,87| of up to 9ae for individual 
Ti-measurements. This strongly limits the accuracy to 
which we can reliably determine ^7^7 1- Another uncer- 
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FIG. 2: Uncertainty in the temperature difference (a) and 
Rb atom numbers (b) introduce uncertainty in the derived 
scattering length. Data for Ti = 10 /iK are shown. The values 
are determined by varying the measured AT and N$7 by small 
amounts dT and dN and performing the fit to the altered 
data as in Fig[T] The influence of the uncertainty dN is small 
compared to that of the uncertainty dT. 

tainty derives from the Rb atom numbers, but this con- 
tribution (5 1 07,87 1 /cW ~ 0.5aB/10 5 ) is small compared 
to the error introduced due to inexact knowledge of the 
Li temperature (FigfJb). We estimate t^T ~ 18 ae and 
£dw ~ laB- The total error in |a7 ; 87| is then 67,37 = 
ye|T+ e^ T + ( 2 dN ~ 19 ob, where % = 5.2 ob is the stan- 
dard deviation of the best fit values averaged over the 
three thermalization measurements. Hence, the result of 
the thermalization experiments is (07^7! = (59 ± 19) ob- 
This result does not agree well with calculations based 
on ab initio scattering potentials [24[, which predict 
0-7.87 = —0.25 ob- With such a value we would not have 
seen any interspecies thermalization on the timescale of 
our experimental cycle. The deviation is likely to be due 
to the very imprecise theoretical knowledge of the scat- 
tering potentials. Probably the good agreement for the 
measured 6 Li- 87 Rb scattering length [J] with the theo- 
retical value presented in [24| is coincidental. 

IV. EVIDENCE FOR COLLAPSE OF THE LI 
COMPONENT 

Upon cooling the Rb to even lower temperatures than 
for the above thermalization measurements, we reach a 
regime where we approach T c for Bose-Einstein condensa- 
tion of Li. Unlike 87 Rb, due to their attractive scattering 
length of aj = — 27.3 <zb, pure condensates of 7 Li are un- 
stable against collapse. However, the kinetic pressure in 
a trapped condensate allows for a certain maximum num- 
ber of condensed atoms [HI, [HI, [25[ . An estimate of the 
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maximum atom number in a stable Li condensate for our 
trapping geometry gives a value of N max « 700. 
Here, dhoi is the length scale of the Li harmonic os- 
cillator ground state. The resolution in detecting atom 
numbers with our absorption imaging is insufficient to 
resolve such small atomic samples. However, observation 
of a sudden drop in the Li atom number for temperatures 
lower than 1.5 /iK (Figj3]) provides a means of revealing 
criticality. Upon reaching the critical temperature for the 
Bose-Einstein transition, Li atoms begin to accumulate 
in the ground state of the trap, forming a small conden- 
sate. The more atoms condense, the stronger the con- 
densate is contracted to a size smaller than the harmonic 
oscillator groundstate [3]. Eventually, the condensate 
becomes mechanically unstable and collapses upon itself, 
thereby strongly increasing the density and thus the rate 
of three body recombinations. Depending on the con- 
ditions, collective decay may also be important or even 
predominant [26l . |27| . All condensed atoms are removed 
from the trap because the energy released in these inelas- 
tic processes is much larger than the trap depth. A new 
condensate is formed on the timescale of thcrmalization. 
With each collapse, this process is repeated, as long as 
the thermal cloud of Li atoms is overcritical with respect 
to Bose-Einstein condensation. This therefore poses a 
permanent loss channel for the much larger number of 
thermal Li atoms. Depending on the thcrmalization rate, 
the overall decay time can be on the order of a hundred 
milliseconds [l8( or many tens of seconds [TB|. In our 
case, we have high trap frequencies and large thcrmaliza- 
tion rates. We therefore expect this process to be rapid. 

While the time-of-flight correction for the Li atoms 
works well for the temperature range used in the ther- 
malization experiments, in the low temperature regime 
we still see a discrepancy between T§7 and T-j. We ob- 
serve T$7 < T-j even though we explicitly took care to 
ensure thermal equilibrium. We hence assume T-j = Tgj 
for this experiment and plot the Li data against Tg7 
(FigJ3|). With the approximate value of the measured 
T c .7 ~ 1.5 /iK we derive a critical atom number of 

N c . 7 = 1.2 • (^f§^) 3 « 6 • 10 5 , which is a factor of 
four larger than the measured atom number. 

If we assume that the losses are really due to a col- 
lapse at T = 1.5 /iK, the small number of measured atoms 
needs an explanation. This number can be experimen- 
tally underestimated if the Li atoms do not remain in the 
cycling transition used for imaging for the whole duration 
of the optical pulse. The cloud therefore becomes trans- 
parent for the imaging beam yielding less optical density 
and lower atom numbers. A complete description of the 
photon absorption process in Li imaging addressing this 
problem in detail is under preparation [28[ . Other groups 
working with 6 Li, which, like 7 Li, suffers from an unre- 
solved hypcrfine structure, image at high magnetic offset 
fields close to the interesting Fcshbach resonance at 822 G 
[29[ | , where the Paschen-Back effect detunes the overlap- 
ping resonances. Full cycling transitions are possible in 
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FIG. 3: The Li atom number drops drastically for tempera- 
tures lower than 1.5 /iK. The resulting critical atom number 
JVc,7 ~ 6 ■ 10 5 is larger by a factor of four than the measured 
value. Possible calibration uncertainties are discussed in the 
text. The dashed lines are to guide the eye. 

this case. 

As an alternative to Bose-Einstein criticality, in- 
creased losses could also arise because of strongly energy- 
dependent three body loss coefficients that derive from a 
Feshbach resonance at low magnetic fields similar to the 
case of Cs |30| . This hypothesis can be safely discarded 
as a mixture of stretched states can in principle not ex- 
perience such interaction resonances. Yet there may still 
be a very small fraction of atoms in the undesired spin 
states such as |2, 1} below our detection limit causing 
residual losses. Even though this cannot be ruled out, it 
is more likely that incomplete repumping during imaging 
causes incorrect atom number determinations and that 
Bose-Einstein criticality is indeed reached [3l[ . 

We have experimentally checked for inelastic losses in 
the mixture in the temperature range between 15 /iK and 
3 /iK. We find that the lifetime of Li in the mixture is 
reduced for colder temperatures, but still remains on a 
timescale of seconds. Going back to the data presented 
in Fig. [31 there, cooling of the mixture from 4 /iK to any 
lower value is done in only 500 ms leading to the observed 
drastic losses. We therefore conclude that inelastic col- 
lisions at thermal densities cannot account for the sud- 
den change in behavior of the data below 1.5 /iK unless 
another loss channel (condensate collapse) opens. We 
attribute this loss channel to be due to the creation and 
subsequent annihilation of Li condensates rapidly deplet- 
ing the observable number of residual thermal atoms on 
a short timescale. 



V. MEAN-FIELD STABILIZATION 

However, the situation may change dramatically in the 
presence of a second degenerate species. The mean-field 
energy can be cither increased or decreased depending 
on the interaction with the second species. If the second 
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species itself forms a stable condensate, as is the case for 
87 Rb, it depends on the sign of the interspecies scattering 
length, whether the collapse of the first species is accel- 
erated or slowed down [32j | . For our system this means 
that for 07,87 < even very small Li condensates are un- 
stable, as for 07,87 > larger meta stable Li condensates 
can be supported. 

To demonstrate this we solve the following coupled 
Gross-Pitacvskii equation for the spherically symmetric 



hid 2 rrik o , 



57,87^;#fc ^fc = M>h 
(3) 

for k,l = 7 or 87. Here g k = 'iirh 2 ak/m,k, 57,87 = 
2-Kh 2 aj 87 m7 +'" 87 and Nk 1 are the atom numbers of both 

*'°' 771777187 ' 

species. 

To find the ground state of the system, we use the 
method of steepest gradient [Hj]. It has already been 
successfully employed to study mixed systems of degen- 
erate bosons and fermions (3J, HH ■ Based on evaluating 
the energy functionals Hk[ipk, tpi], their minimal values 
are found by propagating the wave functions along the 
direction of the gradient -Mf- using small time steps dr: 



4>'k 



tpk + dr 



SHk 



(4) 



Formally, this corresponds to replacing the real time de- 
pendence in the Gross-Pitacvskii equation by an imagi- 
nary time dependence r = it. Under the constraint that 
the normalization of ip' k does not change, which one en- 
sures by explicitly normalizing ipi after each propagation 
step, this method is self-consistent. The final result does 
not depend on the initial trial functions, but a sensi- 
ble choice may significantly reduce the integration time 
needed for convergence. We use the ground state wave 
functions of the respective harmonic trap potentials as 
ipk(T = 0), which are normalized Gaussian functions. 

For simplicity we restrict the calculation to the spher- 
ically symmetric case 14 (r) 



the mean trap frequencies uju 



imfcd^r 2 , where we use 
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The space coordinate was discretized with step size 
dr = 0.1 /Win up to r = 40 yum which corresponds to 
41.9 0^0,87 an d 22.3aho,7- Time was discretized using 
dr = 5 • 10~ 7 s and the propagation was performed for 
6000 time steps. 

For 07,87=0, we get the critical atom number for stable 
Li condensates by increasing N7 to the point where we get 
no numerically stable solution to the above equation. We 
find a value close to the estimated 700 atoms for our trap 
conditions. Setting the mutual interaction to repulsive 
interaction using the measured absolute value to 07,87 = 
+59 ob results in a demixing of the two condensates and 
a substantial increase in the maximum allowed number 
of condensed Li atoms (FigHJ). 

The quoted method only finds the ground state of a 
given system. For situations where collapse occurs, it 



is not suitable to accurately describe the ongoing com- 
plex dynamics. Powerful numerical schemes based on the 
Crank-Nicholson scheme [H, [H, H3] have been developed 
for studying these dynamics, but these go beyond the 
scope of this publication. Here we restrict ourselves to 
the static stabilization of Li condensates in the presence 
of a Rb condensate. Our results may experience mod- 
ifications when considering the cylindrically symmetric 
trap employed in our experiment [38| . Also, different 
gravitational sags for the two species may affect the over- 
all stabilization. This is suggested by experimental re- 
sults presented by D. S. Hall et al. [3{|. They showed 
that the mutual interaction between condensates in the 
two hyperfine states of Rb caused a spherically symmet- 
ric phase-separated state in the case of vanishing rela- 
tive gravitational sag. Introduction of a small gravita- 
tional sag of 0.4 /im, which was smaller than the Thomas- 
Fermi radii of the condensates involved, caused strongly 
damped dynamics and led to a vertically phase-separated 
state with only limited spatial overlap between the two 
condensates. Calculations incorporating the effect of 
gravity will have to be performed to clarify whether or 
not mean-field stabilization then is still significant. In or- 
der to ensure symmetry as much as possible, tightly con- 
fining traps with strongly reduced relative gravitational 
sag should be employed to study mean-field stabilization. 
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FIG. 4: Condensate densities of 1.5 ■ 10 5 Rb |2, 2) (solid line) 
and 6000 mean-field stabilized Li |2,2) atoms (dashed line) 
with trap parameters from our experiment. The interspecies 
mean-field forces the Li condensate into a shell around a large 
Rb nucleus. The resulting density reduction allows for a much 
larger condensate to survive than in the absence of Rb. 



Based on the hypothesis that 07,87 is positive, we per- 
formed a series of experiments where both Li and Rb 
are slowly cooled to about 3 fj,K which is well above the 
experimentally determined value of the critical temper- 
ature for Li condensation. The respective critical tem- 

pcratures behave as T cJ = (jjff) T c , 87 . The 

large mass ratio and the therefore strongly different trap 
frequencies lead to a higher critical temperature for Li 
given our experimental conditions, i.e. T c j w 1.5 /zK 
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while T C] 87 « 0.7 [iK. In order to generate a Rb conden- 
sate without losing too many Li atoms during the evap- 
oration to T c .87, we very quickly (within 200 ms) cool 
the Rb component into degeneracy. The evaporating mi- 
crowave knife remains activated at the end of the evap- 
oration ramp resulting in a slow temperature decrease 
from initially 800 nK down to about 200 nK during the 
900 ms storage time. During the storage stage Rb con- 
densates with atom numbers ranging from some 10 4 up 
to 1.5 • 10 5 atoms are formed. The atom number of the Li 
component drops sharply to less than 10 4 at the end of 
the 900 ms and no bimodal density distribution, the typ- 
ical sign for Bose-Einstein condensation, was observed 
at any time for the Li component. Extraction of the Li 
temperature from the time-of-flight images renders values 
around 1 throughout the duration of the experiment 
while the collisional thermalization model suggests that 
the thermalization is very fast so that AT < 100 nK for 
t > 200 ms. Specifically, the amount of Rb atoms during 
the experiment is sufficient to ensure thermalization even 
with an underestimation of the Li atom number by a fac- 
tor of four. We therefore believe that the Li temperatures 
derived from Li absorption imaging at low temperatures 
are not reflecting the actual temperature of the species. 
Accurate values for very low Li temperatures seem to be 
not achievable under current conditions. A conclusion as 
to the sign of the scattering length is at this stage not 
possible. Two possible scenarios remain: 

(a) The sign is positive, but the relative gravitational 
sag reduces the overlap so that stabilization is sup- 
pressed. The Li condensates stay below the detec- 
tion limit. 

(b) The sign is negative and the Li condensates are 
made even more unstable rendering them unresolv- 
able for our system. 

To reduce the effects of gravity on the overlap of the 
atomic clouds, one should perform these experiments 
again in much tighter confinement. The different trap 



geometry, cigar shaped as opposed to spherically sym- 
metric, has to be considered in the theoretical description 
as well, the get a better quantitative picture of mean-field 
stabilization under our experimental conditions. 



VI. CONCLUSION AND OUTLOOK 

In conclusion, we cooled a mixed gas of two atomic 
species, 7 Li and 87 Rb down to ultracold temperatures. 
We determined the absolute value of the pure triplet in- 
terspecies scattering 1 0,7.87! through measurements of the 
thermalization speed to be 59 ± 19 Bohr radii. The cause 
of slow atom losses experienced by the Li during the cool- 
ing process remains unclear. As Li and Rb atoms in other 
states than |2, 2) are kept below our detection limit direct 
inelastic loss processes can be ruled out. Sudden severe 
Li losses upon cooling below a critical temperature are 
interpreted as the onset of Bosc-Einstcin condensation 
for the Li component. Even though numerical mean- 
field calculations suggest increased stabilization of the Li 
condensate in case of repulsive interspecies interaction, 
Li condensates in the presence of a Rb condensate were 
not detected even in the case of large Rb condensates. 
Such an observation directly determines the sign of the 
interspecies interaction to be positive whereas from the 
current observational viewpoint no decisive conclusion as 
to the sign of 07.87 can be made. The effects of the dif- 
ferent trap geometry and of gravitational sag need to be 
taken into account for more accurate quantitative predic- 
tions. More tightly confining traps that reduce the effect 
of gravity may then prove to be more suitable to study 
mean-field stabilization. A detailed study of the imaging 
process of Li at low magnetic fields is under preparation 
addressing the atom number inconsistencies [28j. 
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